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^sg ' Abstract. We present results from magneto-hydrodynamic jet simulations 

in the magnetically dominated regime and find that they inflate cavities into 
the intracluster medium. Mock Chandra X-ray observations of our simulations 
show a strong resemblance to real cavities observed in X-ray images of clusters, 
both in terms of their morphology and thermodynamic structure. An analysis 
Q . of the evolution of bubble sizes in the multi-cavity system Hydra A, as well as 

Vj ' in a large sample of 64 cavities in 32 clusters shows that bubbles tend to expand 

C/j I much faster than expected in the purely adiabatic regime. Instead, we find that 

C^ . the bubbles follow more closely a trend predicted by our current-carrying jet 

models. 
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1. Introduction 



^^ ■ The cooling time of the hot gas in the center of clusters is significantly shorter 

than the age of the cluster. Thus, the gas should have already cooled completely, 
dropped out of the X-ray phase and formed stars. This so-called cooling flow 

K> ', problem suggests that something is heating the gas, and effectively preventing 

^ ■ it from cooling. Chandra observations of clusters reveal common features in the 

5t , X-ray emission that point toward the central AGN to have a significant impact. 

Images often show deficits in the X-ray emission - generally referred to as cavities 
or bubbles - that are very often nicely filled with extended radio emission (see 



Birzan et al.ll2004l . and references therein). 

So far, theoretical models of AGN feedback have focused on simulating 
purely hydrodynamic, kinematically supported jets. However, these simulations 
have so far failed to reproduce the observed X-ray morphology as depressions 
with enhanced rim emission and the thermodynamic structure of the bubbles, 
which often appear cooler in projection. In fact, bubbles in hydro simulations 
are very quickly shredded apart by Rayleigh- Taylor instabilities on the top and 
Kelvin-Helmholtz instabilities on the sides. Several ideas have been put forward 
to prevent this from happening, with two mechanisms showing the most promise: 
viscosity and magnetic fields. Our work focuses on the second idea, in particular 
the extreme case of magnetically dominated bubbles. 
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Figure 1. Qualitative morphological comparison between an actual Chandra 
observation and a mock Chandra image from our simulation. Left: X-ray 
image of the core of Hydra A; Right: Mock image of a simulated current- 
dominated MHD jet. Note the morphological similarities of depressed X-ray 
surface brightness at the location of the jet lobe, surrounded by enhanced rim 
emission. We emphasize that the rather wide width of the jet structure in the 
center is solely due to limited numerical resolution in this simulation. 



2. Mock Chandra Observations from Magnetically Dominated Jet 
Simulations 



We expan d our earlier work on the first MHD simulat ions of magnetically domi- 
nated jets (JLi et al.ll2006l : lNakamura et al.ll2()06l . l2007l ). Our model only relies on 
the injection of non-force-free magnetic fields in the center of the cluster. These 
fields then expand and launch a magnetic tower. While doing so, the fields first 
self-collimate to form a current-carrying jet structure. When the ambient gas 
pressure starts dropping outside the X-ray core radius, the jet column then ex- 
pands, forming the lobe. During the subsonic lobe inflation process, the lobe 
pushes the hot gas outward, and evacuates a cavity into the cluster atmosphere. 
We then produce mock Chandra X-ray observations with the simulation 
tool MARX. Figure [T]3 shows X-ray snapshots of the evolution of one simulations 
with 200,000 counts. Note the depressed X-ray emission at the the location of 
the lobe, surrounded by enhanced rim emission. The magnetic fields successfully 
stabilize the bubbles against disruption. These key morphological features are 
very close to actual cluster observations, shown for Hydra A in Figure [1^. In 
projection the slow, subsonic expansion of the bubbles lets them actually appear 
cooler than their surroundings. To our knowledge, no purely hydrodynamic jet 
simulation has yet been able to reproduce this observational fact. 
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Figure 2. Left: Black colors show X-ray excess in Chandra's Hydra A image, 
when an elliptical model is subtract ed; circles demon strate the location of 
cavities A-F (Figure reproduced from lWise et al.l l2006. with permission from 
the authors); Right: Bubble size Rb as a function of bubble location r for the 
cavities A-F. The lines show several model predictions. Triple-dot-dashed line: 
purely adiabatic expansion for a F = 5/3; dotted line: same, but for F = 4/3; 
solid line: magnetically dominated, current-carrying jet (our model). Note 
that the evolution of the bubble sizes can be explained by a single current 
running through the system, whereas the adiabatic expansion model would 
requi re subsequently st ronger outbursts. This figure has been reproduced 
from' Diehl eFall (pOOTi ). 



3. X-ray Cavity Sizes: Confronting Theory with Observations 

3.1. Model Predictions: Adiabatic vs. Current Dominated 

Purely hydrodynamic bubble models predict a very simple behavior as they rise 
adiabatically in the intracluster medium (assumed to fohow an isothermal /3- 
model). Thus, - assuming the bubble stays intact - keeping pV^ = constant 
gives a prediction of the bubble size as a function of radius r: 
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(1) 



where Rbfl denotes the fiducial bubble size projected back to the cluster center. 
For our current-carrying MHD simulations we get a very different predic- 
tion. As the lobes are magnetically dominated, the size of the cavity is de- 
termined by the point where the external gas pressure balances the internal 
magnetic pressure (oc B'^). This B-field is mainly generated by the current I^, 
and falls off as Iz/R (Ampres law), thus predicting: 



Rb,l, = Rbfi 1 + {r/r, 
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Figure 3. Left: Scaled bubble size Rb as a function of bubble location r for 
a sample of 64 cavities. The lines are the same as in Figure [^b; Right: Same 
plot, but now symbols are colored according to their black hole masses. Larger 
symbols with lighter c olor indicate large r black hole masses. This figure has 
been reproduced from lDiehl et al.l ()2007l ). 



In this magnetically dominated case, the bubble size Rh^ projected back to the 
center is proportional to the current Iz that is being driven through the system. 

3.2. A Perfect Test Case: The Multi-Cavity System Hydra A 

Unfortunately, we cannot observe the same cavity multiple times as it rises in the 
cluster. However, cavities in the X-ray emission of clusters provide a unique fossil 
long-term record of past AGN activity. The multiple cavity system of Hydra A 
is a great example (Figure [2^), consisting of 3 pairs of bubbles at various radii. 
Assuming that each cavity pair was created by a similar outburst, we can then 
apply our formalism and test the predictions of the models. Figure [2b shows the 
radial evolution of the appropriately scaled bubble sizes (circles) . The red lines 
show the predictions for the adiabatic bubbles (for F = 4/3 and F = 5/3), which 
severely underestimates the bubble sizes at large radii. However, the prediction 
from the magnetically dominated bubbles, denoted by the green line, fits the 
data very well. In order to explain the large bubble sizes at large radii with 
purely hydrodynamic models, one has to invoke successively larger outbursts or 
continuous inflation. 

3.3. Statistical Approach: A Large Cavity Sample 

Figure [3] shows the same plot of scaled bubble size vs. radius for a much larger 
sample of 64 cavities. Note that we observe the same radial trend as for Hydra 
A. Cavities at larger radii are overly large compared to predictions from purely 
adiabatic bubbles. However, the trend nicely follows the general slope expected 
for our magnetically dominated model. It is also worth noting that this result 
is independent of adding viscosity to the adiabatic model, as viscosity cannot 
magically inflate the bubble. The grey lines denote lines of constant current, 
spaced by factors of 4, indicating that the width of the correlation can be repro- 
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duced by a narrow range in currents (about a factor of 30). Figure Et" shows the 
same plot but with symbol sizes indicating the central BH mass. We find that 
more massive BHs drive larger currents and inflate bigger cavities. 

While this r esult is suggestive, one has to keep one caveat in mind, namely 
incompleteness. lEnfilin Sz Heina ( 2002 1) have tackled this problem by analyti- 



cally computing the signal-to-noise ratio of a spherical void in a single-temperature 
/3-model atmosphere. As we do not have an automatic bubble detection tool, 
but rather identify cavities manually by eye, it is difficult to estimate how in- 
completeness really affects our sample. Nevertheless, we can definitely say that 
magnetically dominated cavities will be much easier to identify, as they stay 
intact and expand much faster in general. To put a quantitative constraint on 
the nature of bubbles in cluster atmosphere, we need Monte-Carlo simulations 
of bubbles including incompleteness effec ts, which are curr ently underway and 
will be presented in a forthcoming paper ( Diehl et al.l 120071 ). 



4. Conclusions 

We present results from magnetically dominated, current-carrying jet structures 
in clusters, and show that these jets subsonically inflate bubbles into the intr- 
acluster medium. We flnd that the helical magnetic fleld lines support these 
bubbles and stabilize them against Rayleigh- Taylor and Kelvin-Helmholtz in- 
stabilities. 

Mock Chandra observations show that the inflated bubbles morphologically 
and thermodynamically strongly resemble the cavities found in X-ray observa- 
tions of clusters. In particular, we flnd enhanced rim emission surrounding the 
cavities, as well as an effective decrease in temperature at the bubble location. 

An analysis of bubble sizes in the multiple cavity system Hydra A, as well 
as in a large sample of 64 cavities in 31 clusters favors magnetically dominated 
cavities over purely hydrodynamic bubble s. An analysis of i ncompleteness effects 
will be addressed in an upcoming paper ( Diehl et al.l 120071 ). 



Based on the assumption that our model is correct, we find that the current 
flowing in these systems depends on the central BH mass, offering a potential 
way to constrain BH mass from imaging X-ray observations. 
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